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ABSTRACT 

We compute the decay rate for the process r]c pp using an effective hehcity- 
flipping proton-antiproton-gluon vertex which incorporates nonperturbative chiral 
symmetry breaking effects induced by instantons. We fix the strength of the vertex 
by requiring it to account for the screening of the proton's axial charge observed in 
deep-inelastic scattering, and we estimate the size of the instanton effects by assum- 
ing them to depend linearly on the instanton density. We find that, despite a large 
suppression, the instanton-induced process occurs with a sizable rate comparable 
to the observed one, whereas the process is forbidden in perturbative QCD and not 
understood using standard methods. 



Submitted to: Physics Letters B 



* Address after December 1, 1993: CERN, CH-1211 Geneve 23, Switzerland 







The current theoretical understanding of perturbative QCD has improved at the point 
that it is now possible to single out a few instances where the data seem to disagree with 
perturbative computations. Polarization effects, in particular, seem to display several 
cases where conventional perturbative techniques appear to be inadequate. At high en- 
ergy, where perturbative methods should apply, polarized phenomena are controlled by 
the helicity selection rules which follow from the chiral symmetry of the QCD Lagrangian. 
However, as is well known , this symmetry is partly broken at the quantum level: flavor 
singlet chiral symmetry, which is an exact symmetry of the QCD Lagrangian, is broken 
by the axial anomaly. Besides its standard perturbative consequences (such as the pertur- 
bative mixing of the fermion singlet axial current with gluonic operators 0), this allows 
nonperturbative symmetry-breaking effects, such as those required to solve the U(l) prob- 
lem of QCD This suggests that a related symmetry-breaking mechanism may be 
at the origin of the observed discrepancies with perturbative expectations. Explanation 
of these effects may thus lead to a better quantitative understanding of chiral symmetry 
breaking, a nonperturbative phenomenon in QCD. 

In previous work we have suggested a specific scenario where these effects can be 
understood: we have shown that instanton-induced interactions lead to chirality break- 
ing processes which may be described [0 in terms of a family of effective pseudoscalar 
nucleon-nucleon-n-gluon vertices which flip the nucleon chirality in a way which depends 
on the momentum transfer at the nucleon vertex, but not on the energy of the process. 
Even though the strength of the effective coupling cannot be computed, general arguments 
allow to compute different processes in terms of a single coupling. In particular, assuming 
[3 this coupling to be responsible for the observed discrepancy [|1 between the measured 
axial charge of the nucleon and its quark model value allows to compute the single-spin po- 
larization in elastic proton-proton scattering at high energy and small momentum transfer, 
and predict that it satisfies a scaling law 0. 

Here, we will apply the same physical mechanism to the decay of the r^c- The decay of 
this particle into proton-antiproton poses an outstanding puzzle to perturbative QCD ||^ 
pT[] . This particle is the lightest cc meson; its mass (m^^ = 2980 MeV) excludes decays into 
particles with open charm and all of its decays must proceed through gluon emission and 
are Zweig suppressed. The energy scale for these processes is large enough that they ought 
to be mediated by hard gluons, as required to explain their (observed) Zweig suppression, 
and in agreement with the successful description of analogous charmonium decays. 
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Nevertheless, a perturbative computation of the decay rate for r]c — * pp leads to 
a dramatic disagreement with the data: spin and parity considerations imply that the 
final pp state must be in a total spin S = state, whereas, in the limit of exact chiral 
symmetry, helicity conservation forces the final pp state to have S = 1, hence, the process 
is strictly forbidden. Chiral symmetry breaking effects due to the quark masses can lead 
to a tiny decay rate, of the order of few eV |T^. However, the decay is experimentally 
observed to take place with a width F = (12.1 ± 6.1) keV [jl^. Even introducing effective 
diquark constituents in the nucleon wave function does not help: although these would 
allow the decay without need for chiral symmetry violation (because of the presence of 
spin-1 constituents) the computed rates are typically too small by four orders of magnitude 
[p!0| . An explanation of the large observed r]c pp decay rate seems to require more 
fundamental modifications of the usual perturbative QCD: for example, a nonperturbative 
gluonic coupling of the ijc to the nucleon could provide such new mechanism [|I^ , p!5[1 . 

The effective instanton-induced nucleon-nucleon-n-gluon interaction discussed in 
Ref. incorporates nonperturbative chiral symmetry breaking and could thus provide 
a dynamical model for this kind of process. However, the effect discussed in persists 
at arbitrarily high energies provided the limit of vanishing momentum transfer at the nu- 
cleon vertex is taken; here, instead, the momentum transfer is equal to the mass of the 
r]c- Furthermore, the effective action of Ref. has been derived in the semiclassical ap- 
proximation, which is justified if the strong coupling is large, but it certainly fails in the 
asymptotic as ^ limit; also, in this limit any instanton effect should be exponentially 
suppressed. Here we shall discuss these problems, and estimate the modifications to the 
quantitative consequences of the results of Ref. due to these effects. We will see that 
the rjc mass falls within an intermediate energy range where the strong coupling is small 
enough that a perturbative treatment is allowed, but (essentially because of its slow loga- 
rithmic fall-off) nonperturbative effects are still sizeable and observable in channels where 
conservation laws forbid the bulk of the usual perturbative process. Even though large 
theoretical uncertainties are involved in an estimate of instanton effects, we will be able 
to show that their size is significantly larger than any previously computed perturbative 
contribution and comparable to the experimentally observed rate. 

Let us first review the derivation of the effective nucleon-nucleon-n-gluon Lagrangian 
which we shall use . This is based on the observation Q that instantons may contribute 
to the axial form factor Ga{(1^) of the nucleon, defined by the decomposition 
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where Ga and Gp are respectively the axial and pseudoscalar form factors, u\{p) is a 
nucleon spinor with momentum p, mass rriN and helicity A, and q = p' — p for space- 
like processes and p' + p for time-like ones. That the presence of a classical instanton 
background field may contribute directly to GAiq"^) can be understood by considering the 
time component of Eq. (|l]) in the forward direction, i.e., 

{p,A|Q5|p,A) = 2AG'a(0). (2) 

Eq.(0) shows that Ga(0) is the coefficient of proportionality between the nucleon's total 
axial charge and its helicity: in the presence of instantons, the axial charge receives a 
contribution which is not conserved and corresponds to anomalous violation of chirality. 
An explicit mechanism through which such a contribution may arise has been presented 
in Ref. 0, in a simplified model (QCD with a single massless quark flavor and gauge 
group SU(2), and with the nucleon matrix elements in Eq.(|^) replaced by a quark matrix 
element). 

Whereas this kind of contribution can be computed exactly only in simplifled models, 
such as that of Ref. , its size may be estimated by assuming that it is mainly responsible 
for the fact that the experimentally measured value of Ga{0) is very small, G'a(O) ~ 0, 
whereas the quark model would lead to expect Ga{0) ~ 0.6. Such an assumption is 
supported by the model computation of Ref. This determines the value of the instanton 
contribution to Ga(0). The latter, however, may be viewed as the strength of an effective 
instanton- nucleon-nucleon coupling. Indeed, because of the anomaly equation [Q, Eq. (|l]) 
implies 

lim iG A{q^)qtJiU\'{p')Yl5U\{p) = 



(p, A'l I -2Q+ ^ 2imitlJi'y5i)i ] |p. A), 

\ flavors 



(3) 



where Q equals the number density of instantons minus anti-instantons: 

Nf 

Q = -i^g'tre^'^'F^^F,,. (4) 

Hence, the instanton contribution to the flrst term on the r.h.s. of Eq.(^ provides (in 
the semiclassical approximation) an effective instanton-nucleon-nucleon coupling, which 
we may rewrite in the form 

, A) = -imNGT{Q^)ux'ip'h5Ux{p), (5) 
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where G]^^^{q'^) depends only on q'^ (so that G^®*(0) is a universal coupling) 0. Further- 
more, in the semi classical limit, the instanton can fragment into n gluons (the semi classical 
approximation being allowed if n > ^), thereby leading to an effective nucleon-nucleon 
coupling with the production of n gluons, described by the matrix element [1T6|,|6|| 

n 

{p',X'\l[[Ali{k,)]\p, A) = {2X)zmMGTiQ^)nx'{p'h,ux{p) 

(6) 



X 



n 



16n 



^=1 V 9s 



V^uX' A 1 



kf \ 2 



2 

where K2 denotes the modified Bessel function, Qs is the strong coupling [i.e., as = |^), 
p is the instanton radius, which must eventually be integrated over, a^, and Ui are 
respectively the color*, four- momenta, and Lorentz indices of the n gluons (over which the 
product runs), and fj'^ii,. denotes the 't Hooft symbols i]'^]^^. or fj'^i^. of the instanton or 
anti-instanton (which must also be summed over), respectively. 

In Ref . Eq. (^ was used to derive an effective coupling in the limit of small ki by the 
LSZ procedure; with |G^^*(g^)| ~ |G^'^*(0)| ~ 1 from the requirement that the instanton 
contributions to G'a(O) cancels the quark model one in order to lead to the near- vanishing 
experimental value, and p fixed as the average instanton radius. In such case the strong 
coupling is of order 0;^ ~ 1 and the process is semiclassical with a small number of gluons. 
Here we would like to use the coupling Eq. (||) to describe the production of a gluon state 
which couples to the rjc, so that its decay may proceed through the diagram of Fig.l. To 
this purpose, we must make sure on the following points: first, we should establish that the 
process of Fig. 1 exists, in the sense that semiclassical gluons may couple perturbatively 
to the c quark line; second, we should determine how the size of the instanton effects is 
affected by the the extrapolation from g = to ~ rrirf^ . 

The existence of the process of Fig. 1 relies on the slow fall-off of ctg. At rric the 



value of the strong coupling is known with good accuracy from spectroscopy ||T^, and it 
is given by as(mc) = 0.28. This value is small enough that the coupling of the gluons to 
the charm quark may be treated perturbatively, as is borne out by ample spectroscopic 
evidence. The coupling of gluons emitted from the charmonium decay to the nucleon, 
however, will in general be nonperturbative. Because — ~ 3, instanton-induced processes 



* The indices ai in Eq.(^) run over an SU(2) subgroup of the color gauge group, which must 
be embedded into SU(3). 
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with few gluons are semiclassical at this scale. Hence gluons radiated from the r]c decay 
may couple in principle to the semiclassical vertex of Eq. (||). Because the rjc is C-even 
the leading contribution will be given by the two-gluon process of Fig. 1. Whether this is 
quantitatively significant is the question we address next. 

Quite in general, one would expect that instanton effects are weighted by the Euclidean 
instanton action e ; thus even the moderate decrease in coupling from ~ 1 to ~ ^ 
yields a substantial suppression. In order to estimate this suppression quantitatively, recall 
that the instanton contribution to GAi^) is due to the fact that instantons behave in 
a dielectric way with respect to the axial charge, i.e., in presence of a source carrying 
axial charge they give rise to anomalous creation of axial charge anticorrelated to that 
of the source; thus, the effect will be weighted by the probability of interaction with an 
instanton. We will crudely estimate this probability by assuming it to be linear in the 
instanton density. 

Now, because the effective coupling Eq. (H) is due to anomalous axial charge creation 
in the instanton background, it will only receive contributions from instantons of size 
p < ^, because the anomalous particle creation induced by the instanton occurs at a finite 
distance of order* p. It follows that no particle creation is seen if an instanton of radius 
p is probed at a scale q > j^; consequently instantons with p > pc, with pc ~ -^5 do 
not contribute to the effective interaction Eq. (|D and its cognates Eq. (H). Accordingly, 
the instanton density, on which we assumed G^'^* to depend linearly, will decrease from 
its vacuum value uq (which includes instantons of all sizes) to its value n{pc), computed 
including only small instantons with p < pc. Thus we estimate 

GY\0) no 

The suppression due to the decrease in instanton density can now be estimated using 
Eq. and the explicit expression of the differential instanton density* [|7|, p!9|j20f 

2% 



m^) - ~ (7) 



^ = ^[«^(P )] exp 



(8) 



* For example, it can be proven rigorously [18| that particle creation in an instanton-anti- 
instanton valley background disappears if the separation of the instanton-anti-instanton pair is 
smaller than |p. 
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We use here the vacuum instanton density. In general, the quark-quark chirality flipping 
interaction will carry extra powers of p; these, however, are cancelled by corresponding powers 
of p in the quark propagators in the instanton background when computing a matrix element, 
such as that of Eq. @, as it is clear on dimensional grounds. 
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where all the p-independent quantities have been lumped in the constant C. The differ- 
ential density Eq.(P) should then be integrated over all instanton radii. The density of 
instantons with radii up to pc will be given integrating in the range < p < Pc', because 
in the case of interest pc is in the perturbative region the integration can be performed 
by using the perturbative expression for as- The vacuum instanton density in the de- 
nominator of Eq. instead, of course diverges for large instanton radii; this problem is 
cured phenomenologically |T^ by assuming the growth of the instanton radius to be cut 
off at some scale po ~ 1 Fm by a (poorly known) instanton repulsion mechanism. The 
corresponding value of the average instanton radius (obtained averaging with the measure 
Eq. (H) up to po) is of order of (p) ~ ^ Fm. 

In order to minimize the model dependence of the computation of N{q), Eq. (|^), 
rather than using a phenomenological value for no we compute the denominator of Eq. (^ 
by integrating Eq. (|^) up to a cutoff value po ~ 1 Fm. This indeed minimizes the model- 
dependence of the result because then the value of N{q) is controlled by the exponential 
dependence of n{p) on p around p ~ "^^^^5 while it is essentially independent of pre- 
exponential factors. The former in turn is controlled by the value of asirrij^J, which is 
known rather accurately. On the contrary, as we will shortly see, the result is rather 
insensitive to both the upper limit of integration po, and the precise value of as{po^), 
which is large, as{po^) > 1- 

Due to sensitivity of the results to the perturbative running coupling, we use the 
next-to-leading form 

d-rr r ft. In In ^1 

(9) 



l3o\nt 



/3ilnlnt 
/3^1nt 



where [3q = 11 — 2n//3, [3i = 102 — 387ij/3 and in the perturbative regime t = ln((5^/A^). 
In order to compute the suppression Eq. (0), we must introduce an infrared interpolation 
for the strong coupling a^, so that at large the perturbative behavior is reproduced, 
while at small the strong coupling ctg saturates to a constant value. To this purpose 
we set t = 5 + ln((5^/A^). For the perturbative QCD scale we take the value (with n/ = 4) 
A = 263 [|2l|, which yields cts(mc) = 0.31, in good agreement with spectroscopy, whereas 
5 is fixed by the value of as in the infrared. For example, imposing Il|] a^(lFm"^) = 1 



we get 6 = 3, while larger infrared values of as can be obtained by reducing the value of 
5 (which must anyway satisfy 5 > 1 to insure saturation). 
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With this form of the strong couphng*, 5 = 3, and taking the cutoff radius of the 
vacuum instanton density to be po = 1 Fm, we get N'^irrir^^) ~ 3.9 x 10~^. Varying the 
infrared value of by an order of magnitude as weU as the value of po from 1 Fm to 10 Fm 
this determination varies at most by a factor two, thus displaying the infrared stability of 
our estimate. Even though there is (as expected) a substantial suppression of the coupling 
Eq. (H) due to the decrease in instanton density, the strength of the coupling is still large 
enough to lead to a sizable decay rate, as we show next. 

We proceed therefore to the computation of the diagram of Fig. 1. The coupling of 
gluons to the charm quark line is assumed to be given by perturbative QCD; the charm 
quark-antiquark pair then hadronizes with a wave function that we assume to have the 



static (nonrelativistic) form, p3|,p!0 



4v^ 



(10) 



where i,j are colour indices and ± denote the quark helicities. The value of the radial 
charmonium wave function in the origin, -R(O), can be fixed computing the width of 
the electromagnetic process r]c 



r(^c ^ 77) 



77 and comparing to the experimental value |]13 
64 a2 



27 ml 



■|i?(0)ni-3.4^) = (6.6i2-t) keV, 



-2.1) 



'ID 



which yields i?2(o) = (0.67tO-J4) (GeV)^. 

The amplitude for the decay process is then given by 



A 



A„A, = 7ri?(0) Ma,,a,;++(A; = 0) - Ma,,a,; 



.(fc = 0) 



where M are the elementary helicity amplitudes for the process cc pp: 

d^ki d'^k2 



M 



(27r)4 (27r)4 

RZ xM k) T':i^^^{k^M.p) {2'nf5^-^\ki + ^2 - 2c), 



(12) 



(13) 



The computation of the integral over instanton radii has been performed by taking the 



appropriate number of flavors at any scale, i.e., n/ = 3 for i < rric, rif = A for rric < - < 



rub 



and so forth, and suitably updating the value of A [22|. The result turns out to be essentially the 
same as that obtained setting = 4, due to the rapid falloff of the instanton density at large p, 
as well as the insensitivity of the value of N{q) Eq. (^) to the precise form of the coupling in the 
infrared. 
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and we have set c = (mc, 0), rric = and k = {ki — A;2)/2, while R^^ and T^'^ give 

respectively the perturbative part of the amplitude, and the nonperturbative coupling to 
the proton. Explicitly, the perturbative contribution is given by 

KlxM^ k) = --A.(c)7'^^|^7^«A.(c), (14) 

where we have Wick-rotated to Euclidean space where the instanton coupling is defined. 
The nonperturbative coupling is constructed by taking n = 2 in Eq. and calculating 
the color trace, which, due to the form of the wave function Eq. ( |ToD amounts basically to 
projecting out the color singlet component, with the result 

(15) 



X 



ki ■ k2g^'' - ^{k'^k!^ + k^ki^) + (2Ap)e-'^^^A;i«A;2/3 



In Eq. (^5]) we have set p' = ki + k2 — p, and 

1 



$(fc)=4(l-^i^2(p|A:|)pV). (16) 



Using Eq.s (|I4D and ( [TSD in Eq. (|T3D we get, through a somewhat tedious but straight- 
forward computation 

Mx^, Xp A„ A, = ^^^GA{m^J 4:ml {4X^Xp) Sx^x.Sx^Xp 



(17) 

(27r)4 A;4(2c- A;)2[(c- A;)2 + m2 
The loop integration in Eq. (|T^) converges, as one would expect due to the effective nature 



of the interaction (H), because the function $(/c), Eq. (|T6|) , behaves as 



2l2 



(j){k) ^ p k 



(18) 



thereby cutting off the ultraviolet tail of the integration while being finite in the infrared. 
Using the result Eq. ([T^) in the expression (|T^ of the decay amplitude and supplementing 
the required kinematical factors the width for the process is found to be 

nvc^pp) = Iq^2 E i^a„a,i'. (19) 
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In order to determine the value of F, Eq. (p!9[) , we must still fix the values of the instan- 
ton radius and of the coupling constant in the nonperturbative portion of the amplitude, 
Eq. (|T3p. In principle, all p-dependent quantities ought to be included in the averaging 
over instanton radii which also provides the decrease of G]^^^ according to Eq. (^; this, 
however, would unnecessarily (to the level of accuracy of this estimate) complicate the 
evaluation of the loop integral in Eq. (|T^. Hence, we evaluated instead the integral by 
fixing p = Using the value quoted above of N{mnJ we get thus finally 

r ~ 0.8 keV (20) 



We can get a more direct handle on the ingredients which enter the estimate Eq. ( poD 
by noting that, because of the behavior Eq. ( |T8D of the function $(/c), the loop integration 
is effectively cut off in the ultraviolet at momenta of order ~ ^. If accordingly we 



neglect the contribution to the integral from momenta \k\ > -i, using the low- momentum 



form of Eq. (^8|), the integration may be performed analytically, with the result 



Kx\I\^x\N\^x4S^. (21) 



where all kinematical factors have been lumped in 



K = |i^(0)|^ (22) 



while / is the result of the loop integral. 



CI 

1 



"^^J "^'^ (2c - A;)2 [(c - A;)2 + ml] ^^3) 

= P TT In r-^ 

and N provides the overall suppression due to the instanton density, Eq.(^. Numerically, 
K ~ 0.3 (GeV)^ while (with p'^ = m^J / ^ 0.2 (GeV)"'^ and A^^ ^ 3.9 x lO'^, which 
indeed reproduces the result Eq. (pO]). 

It thus appears that the decay rate induced by the chirality-flipping interaction is 
kinematically of order of several hundred MeV, but it is then reduced due to the large 
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suppression of instanton effects to the order of the keV, in good agreement with the exper- 
imental value. Clearly, this is a crude estimate of the order of magnitude of T{r]c pp) 
induced by instantons. The main uncertainties involved in this estimate enter in the de- 
termination of the suppression factor N{q) Eq. (0), mostly in the exact dependence on the 
cutoff instanton radius and in the precise value of the strong coupling. For example, if the 
value Pc = f ^ were used (appropriate for the instanton-anti-instanton valley of Ref. [1T8[|), 
rather than Pc = then the value of N(q) would increase by an order of magnitude. Anal- 
ogously, varying A between 200 and 300 MeV the value of N{q) fluctuates by about one 
order of magnitude. These uncertainties also affect the effective value of p which should be 
used in the computation of the loop integral Eq. (^7D , on which the loop integral depends 
as / oc /9^, according to Eq. (pSf). As discussed above, a smaller uncertainty comes from 
the infrared behavior of the instanton density. Finally, a minor uncertainty enters in the 
experimental determination of R{0). In short, most of the uncertainty comes from the 
perturbative behavior close to the mass of the rjc, which is under theoretical control and 
may be reduced by a more accurate treatment and a better determination of A. The result 
turns out to be in remarkable qualitative agreement with experiment. 

The instanton-induced coupling discussed here leads also to decay of the first radial 
excitation of the ry^ according to the same mechanism, i.e., to the process ry^ pp. This 
has not been observed yet, even though data are expected to become available soon. If 
these processes proceed mainly through the instanton-induced mechanism discussed here, 
we expect the value of r(r7^ pp) to be smaller than that for the rjc. Indeed, even 
though the larger value of the 77^ mass [rrirj'^ ~ 3.6 MeV) leads to an increase of the phase 
space (compare Eq. ([2^)), this is largely compensated by the very rapid decrease of the 
instanton suppression factor Eq. (^ (and to a lesser extent, the loop integral Eq. (|3l)) as 
the characteristic scale of the process increases. Specifically, we get 

g_ r(r?-^PP) ^ \Rmv'^)\' 

-Tiv^^pp) - ■ |i?(0)(r?,)|2 ^ ' 

The precise value of S depends on the ratio of the radial wave functions in the origin. 
Assuming for instance this ratio to behave as that of successive s-wave radial excitations 
in a hydrogen-like potential, then i?2o(0)/-Rio(0) = l/"\/8, leading to ~ 3 x 10~^. 

In sum, we have estimated the magnitude for the width of the decay r]c — * pp using 
an instanton-induced effective Lagrangian to model the chiral symmetry breaking which is 
necessary in order for the decay to take place. The decay proceeds through the emission by 
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this effective Lagrangian of a gluon pair which is soft enough to be treated semiclassicaUy 
at the proton vertex, but hard enough to couple perturbatively to a charmed quark hne. 
The result turns out to be remarkably of the correct order of magnitude, whereas all 
previously known mechanism to describe exclusive processes within QCD had failed (by 
several orders of magnitude) to give the correct result. The main uncertainties involved 
in this estimate come from the determination of the strength of the instanton effects at 
the scale considered here, where, even though nonperturbative effects cannot be neglected, 
one is already well within the perturbative region. Since these effects can be brought 
under theoretical control, this provides a strong motivation for a systematic investigation 
of instanton-induced helicity effects at medium-high energy of which Ref. (for forward 
elastic scattering) and the present work are the first steps. 
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FIGURE CAPTION 

The diagram which leads to the decay r^c PP- The blob indicates the effective 
instanton-induced interaction of Eq. (H) . 
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